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SYMBOLS AND CONVENTIONS 


aspect ratio = 2 

mean aerodynamic chord (ft.) 

teil mean aerodynamic chord (ft.) 
elevator mean aerodynamic chord (ft.) 
center of gravity 


lift coeff. = & 


q8 
dor 
slope of lift curve = — +4 
dx 
d 
lift coeff. of horizontal tail due to elevator deflection - (τε ) 
t 


Stick fixed static stability = E 


downwash lag derivative - dC 


d (Do) 


dC 
in pitch derivative = m 
damping p 17555 


elevator power derivative = ACm 
ds 





differential operator with respect to dimensionless time 
density altitude (ft.) 


2 
non-dimensional inertia about y axis = 2 (4) 


A- NC 
moment of inertia of airplane about y axis (slug ft.^) 
radius of gyration K, 2. ix (ft.) 
tail length (ft.) (dist. from airplane c.g. to a.c. of tail) 


mass of airplane = Ξ (slugs) 


mean aerodynamic chord (ft.) 





Το 


Νο stick fixed neutral point (% M.A.C.) 

Na normal acceleration - positive up 

g dynamic pressure = dev” (lb. per sq. ft.) 
S wing area (sq. ft.) 

S4 tail area (sq. ft.) 

Se elevator area (sq. ft.) 

T.R. taper ratio 

V velocity (m.p.h.) 

W gross weight of airplane (lbs.) 


Xc4. location of c.g. ($ M.A.C.) 





GREEK TERMS: 


A 


δ 


B, 


angle of attack of wing (degrees) 

elevator angle (degrees) 

downwash angle (degrees) 

tail efficiency factor 

density of air (slugs per cubic ft.) 

rate of change of downwash with angle of attack 
= TEE h 

aerodynamic time psc (sec.) 

airplane density factor = mM. 

ϱ 86 

phase angle, minus for lag angles (degrees) 


angular frequency (rad./ ec.) ; or (a, e) if specified 


FLIGHT TEST DATA SYMBOLS (Fig. 12-23): 


ο 
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Δ 
n 
o 
Y 


Flight Test 71 
Flight Test #2 
Flight Test #3 
Flight Test fl 
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Flight Test # 





AIRCRAFT SPECIFICATIONS 


Ryan Navion N5113K 


Average Weight as flown 
Wing Area 

Horizontal Tail Area 
Elevator Area 

Aspect Ratio of Wing 
Aspect Ratio of Tail 
Wing Span 

Tail Span 

Wing MAC 

Tail MAC 

Elevator MAC 

Dist. from c.g» to Tail a.c. 
Taper Ratio - Wing 


Airplane Inertia, Y axis 


2825 lbs. 
181.2 8Qe 


13 Sq. 
15.0, 89 
6.04 
“4.02 
33.38 ft. 
13.17 ft. 
Saf) cae 
3.31 ft. 
1.20 tt. 
15.0] ft. 


1.83/1 


ft, 
fto 


fto 


3,000 slug ft.? 
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SUMMARY 


The purpose of this investigation was to determine the feasi- 
bility of obtaining reliable longitudinal frequency response data for 
an airplane using simplified techniques, equipment and instrumentation. 
The particular method involved sinusoidal oscillation of the elevator 
at various frequencies and measurement of the airplane's response. 

The instrumentation was designed and built with the basic 
premise of simplicity. The elevator was oscillated mechanically, and 
the response was recorded by a small, lightweight SFIM recording oscíl- 
lograph. The tests were conducted at two velocities with three differ- 
ent airplane center of gravity locations and the data reduced to general 
frequency response form. For one of the test cases a comparison was 
made between the theoretical and experimental responses an analytical 
expression for the transfer function based on the experimental data was 
also determined for this case. 

This approach proved entirely feasible for obtaining frequency 
response of the Ryan Navion as repeatable results were obtained with 
little random.error. The same instrumentation scheme could likely be 
employed in many current production aircraft as the instrumentation 
space and power requirement 1s small, This type of testing is however 
restricted to flight vehicles in which steady state conditions can be 
achievedo 

If stability derivatives are to be obtained from the results 
of this type investigation, a limitation on the degree of simplicity 
may be imposed. 
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INTRODUCTION 


Dynamic flight testing of airplanes has until recently been 
concerned almost completely with the qualitative approach of determin- 
ing whether the airplane was dynamically stable or unstable. With the 
growing development in the field of automatically controlled and stabi- 
lized aircraft, a quantitative determination of dynamic response char- 
acteristics becomes a greater concern. This information is required 
for the control system designer so that he can consider the airplane 
with its aerodynamic transfer function as another link in the overall 
system. 

The determination of this transfer function or frequency re- 
sponse for an airplane may be accomplished in a number of ways. First, 
it could of course be derived analytically, but this introduces all the 
attendant inaccuracies and difficulties associated with the equations 
of motion and stability derivatives. A more direct approach is to de- 
termine frequency response from either transient or steady state dynamic 
flight testing. The transient method has the disadvantage of requiring 
complicated recording instrumentation, and the resulting data is quite 
difficult to analyze. Its principal advantages result from the economy 
of flight time requíred to obtain data, and from the lack of necessity 
for steady state flight conditions. Steady state dynamic flight test- 
ing results in very simple data analysis with the possibility of less 
complicated instrumentation, It does however require considerably more 
testing time and steady State conditions are necessary; these demands 
can limit the applicability of the method. 


Longitudinal frequency response data may be obtained for an 





airplane through steady state dynamic flight testing utilizing a forc- 
ing function upon the airplane's elevator. This forcing function con- 
sists of a sinusoidal oscillating elevator input of controllable ampli- 
tude and frequency. The frequency response data may then be obtained 
by recording the airplane's response in normal acceleration, pitch rate, 
etc. for a steady state condition. The presentation of the frequency 
response data is accomplished by measuring and plotting vs. frequency 
the amplitude ratio of that particular function of airplane response 
desired to the elevator input and the phase angle between them. 

Tests of this sort have been conducted by the Cornell Aero- 
nautical Laboratory employing a system utilizing an autopilot with a 
sine wave generator to give the purest possible sine wave oscillation 
to the elevator. This led to considerable complication in the instru- 
mentation which made this particular approach impractical for other 
than scientific types of testing. 

This report was conducted to investigate the feasibility of 
determining frequency response of airplanes by the same type steady 
state dynamic testing used in the Cornell Aeronautical Laboratory tests. 
However more simplified instrumentation, applicable to production type 
airplanes, was used. The system herein discussed involved a mechanical 
elevator oscillation scheme and the SFIM recorder. The flight tests 
were conducted in a Ryan Navion at the Forrestal Research Center, 
Princeton University. 





13ο 


INSTRUMENTATION 


General. 

Investigation of the frequency response of an aircraft poses 
two instrumentation problems. First, the aircraft response must be 
measured and recorded and, second, a sinusoidal elevator excitation 
must be supplied. In addition, for this thesis simple instrumentation 
was a primary consideration. The SFIM A-20 flight recorder is an ex- 
tremely small recording oscillograph that met the above requirements 
and was selected for use in this investigation. This instrument is a 
2l volt D.C. machine, so it was decided that only compatible 2|, volt 
D.C. components would be employed. 

The oscillating device was built around a 2l volt surplus 
antenna reel motor. The choice of where to mount the motor and how to 
convert the rotary motion to an accurate sinusoidal motion was also 
dictated by simplicity. A drive shaft and eccentric arm appeared to 
be the most simple and accurate conversion. Mounting the device in the 
tail and connecting it directly to the elevator was most desirable but 
proved impossible because of space limitations, inaccessibility, and 
Safety of flight considerations, Finally, it was decided to mount the 
unit immediately aft of the front seats with the driving arm extending 
between the pilot and copilot and connected directly to the base of the 


control column. This has proven a most satisfactory arrangement, 


Components. 
The Aircraft. (Fig. 1 


The aircraft used in this investigation was a standard 


Ryan Navion number N5113K, The only modification to the air- 
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plane was the addition of a sliding weight mounted along the 
longitudinal fuselage axis which was used to control the air- 
oraft's center of gravity position. This is described in 
detail in Princeton Report 275. 

Electrical Power Supply. | 

A 2l volt D.C. two wire system was used. The aircraft's 
12 volt battery was connected in series with a second 12 volt 
heavy duty battery to provide a 2l volt source of power. 

This power source supplied a bus bar barrier strip mounted on 
the starboard deck aft of the right rear seat, 
SFIM A-20 Flight Recorder, (Fig. 3) 

This is a very small (6.5" x 1.75" x hl") ten channel re- 
cording oscillograph. Recording of various traces is made on 
a photographic paper strip by beams of light supplied by the 
filament of a lamp located in the recorder. This image is 
projected by mirrors of the measuring instruments which are 
installed inside the recorder. The recorder incorporates a 
timing unit which is regulated by a clockwork escapement. 
Four sockets mounted two on each lateral flange are provided 
to hold measuring instruments. In this investigation only 
four channels were utilized, and they recorded elevator de- 
flection, normal acceleration, pitch rate, and time. 

The SFIM A-20 flight recorder proved to be a most excel- 
lent piece of equipment. The only difficulty experienced in 
its use was that the film speed of the machine is too slow, 
approximately 1 ft./min., and this made some high frequency 


traces difficult to analyze for phase angles. The recorder 
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holds twenty feet of film which proved to be an adequate 
amount for each flight. 
Normal Accelerometer. 

The accelerometer used was a SFIM J531 instrument mounted 
internally as channel 2 in the recorder, The detecting ele- 
ment is a suspended mass. Range is 0 to 2 g. It is air 
damped and has a natural frequency of 25 cycles per second. 
This is well above the range of investigation so no instrument 
resonance was encountered. 

Rate of Turn Gyro Transmitter. (Fig. h) 

The instrument used was a SFIM I 131 rate of turn gyro 
transmitter. This device transmitted electrically measure- 
ments of angular speeds, The measuring element is a 2l volt 
D.C. gyroscope. The transmitting element is a D.C. potenti- 
ometer. Speed of rotation under 27 volts is 9000 RPM, 
Natural frequency of this component is 5 cycles per second. 
The range of this investigation is between 0.2 and 0.5 
cycles per second, so the dynamic characteristics of this 
instrument should be accounted for. Range of the instrument 
is 69/sec. The instrument is & small cubic box measuring 
3.25 inches side length. 

The signal transmitted by the potentiometer is received 
and inscribed on the recording film by a ferromagnetic rati- 
ometer type E 521 mounted internally in the recorder as 
channel 1. 

Elevator Deflection, 

The elevator position was measured by a Helipot contim- 

ous rotation potentiometer which was located in the tail of 





the aircraft. <A rigid arm was mounted to the hinge line of 
the elevator and a string-spring combination used to transmit 
elevator motion to the potentiometer shaft, The potentiome- 
ter was powered by the 2l volt D.C. system and the output 
signal was received and inscribed on the recording film by & 
ferromagnetic ratiometer type E 521 mounted internally in the 
recorder as channel lj. 
Oscillating Device. (Figo 

The elevator oscillating device is a very critical item 
in an investigation of this type. In order to make the | 
flight results of value, the device must drive the elevator 
sinusoidally at & given frequency and amplitudeo Theoretical 
Studies of the airplane's modes of motion indicated that the 
proper range of frequencies to be investigated was 12 to 150 
cycles per minute, Examination of the elevator hinge moments 
and gearing showed that less than 0.1 HP would be required to 
deflect the elevator 10 degrees at a "q" of 50 p.s.f. The 
motor that most nearly met the above requirements at a moder- 
ate cost was a surplus 2l volt D.C. antenna reel motor of 
1/8 H.P. @ 1000 RPM geared to 150 RPM output. This motor 
has two field windings for forward and reverse operation, 
However, in this installation the motor was required to run 
in only one direction and resourceful utilization of the 
second field winding provided a built in tachometer generator, 
The current induced in the second field winding by the arma- 
ture was measured on an A.C. voltmeter which was used as a 


tachometer. The motor speed is controlled by a variable 
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resistanoe in series with the motor. This resistance con- 
sists of a five ohm rheostat and a 3-1/3 ohm resistor which 
may be added to the circuit by a switch. 

The motor is mounted with the shaft parallel to the 
airplane's lateral axis and drives a small flywheel through 
a universal joint. The outboard face of the flywheel is 
drilled and tapped at four radii which represent elevator 
deflections of approximately 2.59, 59, 7.59, and 10%, Eleva- 
tor deflection is achieved by mounting a 5/8" grooved pin at 
the desired radius in the flywheel face. A similar pin is 
fixed to the base of the control colum. A drive shaft con- 
nects the flywheel pin and the control column pin and con- 
verts the circular motion of the flywheel to the sinusoidal 
motion at the control column. The ratio of flywheel radius 
to drive shaft length is less than 1/50 at the maximum eleva- 
tor deflection. The drive shaft has collars with spring 2 
clips on both ends. The spring clips fit into tbe groove in 
the pins on the flywheel and control column. The pins are 
made of mild steel and the collars are made of bronze, hence 
& proper bearing surface is achieved, In addition the drive 
shaft incorporates a turnbuckle arrangement for trimming the 
aircraft and sleeve and setscrew arrangement for gross 
length adjustment. 

The oscillating device is mounted on a l/h" boiler 
plate base which is attached to the sliding weight housing 


immediately aft of the front seats. 
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Instrument Control Panel. (Fig. 2) 
An instrument control panel was mounted on the forward 


face of the rear seat instrument console. This unit contains 
an A.C. voltmeter which is used as a tachometer, a D.C. volt- 
meter for measuring battery potential, the motor rheostat and 
step resistance, and on-off switches for 2l| volt power, SFIM 
recorder, rate gyro power and signal, and elevator deflection 
Signal. A fast/slow film speed switch was placed in the 
panel to accommodate both film speeds available in the SFIM 
recorder. Since only the fast speed of the present recorder 
was utilized, this switch was left unconnected. In addition, 
two four-prong receptacles are provided for electrical con- 
nections to the SFIM recorder and motor, A three-prong re- 
ceptacle is provided for electrical connections to the rate 
gyro. The instrument panel is powered by connections to the 


2); volt D.C. bus bar barrier strip which is described above, 


Calibration. 

Four items required calibration in this instrumentation sys- 
tem. They were the rate gyro, the normal accelerometer, the elevator 
deflection indicating system, and the tachometer. The aircraft air- 
speed indicating system was not calibrated because all previous calibra- 
tions showed that I.A.5. could be used as C.A.5. with less than 2.0 MPH 
error for normal aircraft configurations. 

The rate of turn gyro was calibrated on a pendulum which con- 
sisted of a 3/1” diameter aluminum shaft and a 12" diameter instrument 


pano The pendulum was pinned to the overhead truss structure in the 
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hangar. See Fig. 6. Overall length of the pendulum was 120.75". The 
rate of turn gyro and the SFIM recorder were mounted on the instrument 
pan and power was supplied by a 2l volt battery cart through wires 
which extended down the shaft from the hinge pin to the instruments, 
Theoretical calculations, based on a pendulum length of 120.75", gave 

a period of 3.505 seconds. The actual period as determined from the 
SFIM recorder timer trace was 3.6 seconds. This discrepancy was the 
result of theoretical calculation based on the overall pendulum length 
rather than the length to the c.g. of the system. Turn rate was calcu- 
lated from the relation 6 = Ow cos wt so that 5 = Bw 
where w was actual pendulum frequency and 9 was measured angular 
displacement from zero. Angular displacement was varied in five steps 
from .083 radians to .1823 radians. The pendulum proved to be isochro- 
nous for these amplitudes. Pitch rate gradient was found to be linear 
and equal to 0.1222 radians/sec. per inch trace amplitude cycle. 

The accelerometer was calibrated by securing the SFIM re- 
corder to an adjustable tilt head band saw table. The table angle was 
measured with a propeller protractor and a record was taken in steps 
from -709 to / 709 tilt. Normal acceleration was calibrated as the 
cosine of the tilt angle with the zero position equal to 1.0 g. Power 
for the recorder was supplied by a 2l volt battery cart. Normal accel- 
eration gradient was found to be linear and equal to 0.68 go per inch 
trace amplitude cycle. 

The elevator deflection was calibrated by measuring elevator 
deflections between -10° and 410° with a propeller protractor and re- 
cording these deflections in one degree steps. For this calibration 
the SFIM recorder was installed normally in the aircraft instrumenta- 
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tion system and powered by the 2l volt D.C. power source. The elevator 
deflection gradient was found to be non-linear at the extreme deflec- 
tions, but was linear within a 5? deflection range about the required 
elevator trim position. The gradient was found to be 5.97 degrees per 
inch trace amplitude cycle. 

The tachometer was calibrated in flight by timing the motor 
RPM at various resistance settings and then assigning RPM values to 
measured induced voltage readings on the A.C. voltmeter. Accuracy of 
the tachometer system allowed only order of magnitude readings, but 
this reading proved to be of great value ín the flight test program. 
The actual frequencies were determined later from the SFIM recorder 


timer record. 


Instrumentation Deficiencies. 

Experience uncovered two rather serious deficiencies in the 
instrumentation system. The first of these was the speed of film 
travel in the SFIM recorder. Two speeds are available in this model, 

l mm/sec and 5 mm/sec. Although the faster speed was employed through- 
out the investigation, considerable difficulty was experienced in ob- 
taining phase angles at high frequencies. Low frequency phase angles 
could be found with some accuracy, but it is felt that high frequency 
results are unreliable due to difficulty in data reduction. This film 
speed defect created no amplitude measurement problems however, and 
amplitude response is considered satisfactory. Unfortunately, there 
was no method of increasing film speed in this recorder and this short- 
coming remainded uncorrected. 


The second deficiency was in the oscillating device. Unfort- 
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unately, motor torque varied directly with voltage, and while the 
motor could be made to run slowly, it produced too little torque to 
operate the oscillating system. This inadequacy was eliminated by 
the addition of a 4.5:1 reduction gear and a second drive shaft. Two 
mounting positions for the motor were established so that the oscil- 
lating device could be driven by either the low or high speed shaft. 
The shaft connection can be changed in approximately 15 minutes. No 
further difficulty was experienced with this problem. 

A third, rather disconcerting defect was discovered on the 
initial shakedown flight when the oscillating device was being tested. 
The spring clip releases on both front and rear collars of the drive 
shaft linkage were designed as quick release devices and had been 
satisfactorily tested on the ground. They did not work in the air 
however, and the airplane, instrumentation, and authors were very 
nearly destroyed before manual control could be restored. After this 
experience a T handle set screw arrangement was d SEDO ESA in the 
drive shaft linkage. See Fig. 5. This proved to be of aid in start- 
ing the oscillations in addition to giving peace of mind as an overall 


safety device. 
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FLIGHT TEST PROGRAM 


The flight test program was designed to obtain data for síx 
flight conditions at 5000 feet density altitude. These flight condi- 
tions were for two airspeeds, 100 MPH and 110 MPH, at three center of 
gravity positions; 25.0, 29.85, and 34.60% ma.c. As experience was 
gained in data recording procedure, the 20 minute film duration of the 
SFIM recorder proved to be sufficient time to investigate all six of 
these conditions on each flight. Due to the change in motor installa- 
tion required to investigate the entire frequency range, it was de- 
cided that only one frequency range should be investigated on any 
flight. Low frequency data was obtained first and after sufficient 
repeatable results were obtained the motor waS remounted and the test 
program repeated to obtain high frequency data. 

Each flight test crew consisted of a pilot and a copilot- 
recorder, and the following procedure was adhered to on each flight. 
Prior to take off the power source was checked for a 25 volt minimm 
indication and the instrumentation system was tested. While climbing 
to 5000 ft. density altitude the rate of turn gyro was uncaged, the 
sliding weight was shifted to obtain the desired c.g. location, the 
instrument check list completed, and the drive shaft linkage connected 
with the set screw unlocked so that the pilot retained normal control 
of the aircraft. When at altitude the pilot adjusted power and eleva- 
tor trim tab for steady level flight in the desired flight condition, 
With the drive linkage still unlocked the copilot turned on and ad- 
Justed the speed of the oscillating device, The pilot then locked the 


set screw on the drive shaft linkage thereby completing the connection 
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from the device to the elevator. It was necessary for the pilot to 
lock on in a position that caused the elevator to be oscillated about 
the elevator trim deflection, otherwise a somewhat uncomfortable nose- 
up or nose-down trend accompanied the normal oscillation. A satisfac- 
tory lock on usually required several attempts. After the gross ad- 
justment had been achieved by a satisfactory lock on, airspeed and 
altitude were observed and final adjustments accomplished by the turn- 
buckle linkage trim systemo Patience was required in this process but 
skill was acquired with experience. After all adjustments were com- 
pleted, the copilot turned on the recorder and at least ten steady 
State cycles were recorded. The recorder was then turned off and the 
linkage unlocked. The copilot made an entry in the data log of air- 
speed, altitude, tachomster reading, and any additional pertinent com- 
ments. He then set the rheostat for a different frequency and the 
lock on and recording procedure was repeated until all desired fre- 
quencies in that condition were investigated. A new condition was 
then effected by changing the c.g. location or airspeed. The most 
efficient method of covering all conditions was determined to be to 
start the runs at 100 MPH in the aft c.g. condition, then go to 140 
MPH at the same c.g. location, next to shift to mid c.g. while main- 
taining 140 MPH, etc. with the final run made at 10 MPH with c.g. 
forward. At the completion of the tests, the rate tyro was caged and 
the instruments secured. 

Flights 1 through were investigations in the low frequency 
range» After flight the motor was remounted and flights 5 and 6 
were investigations in the high frequency range. On flight 1, a 5/8" 


foam rubber pad was placed under the SFIM recorder and the rate gyro, 
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This was done to eliminate trace distortion due to airframe vibrations, 
On flight 2 the sponge was removed and both flight records were studied. 
There appeared to be no vibration distortion in the flight 2 record so 
the sponge pad was not used on subsequent flights. In addition, it is 
felt that some inaccuracy may result in the normal acceleration record 
due to the cushioning effect of the foam rubber pad if it is used. 

Rough air makes this type testing prohibitively difficult. 
Moderate turbulence was found to continuously excite the phugoid mode, 
as evidenced by slow altitude and airspeed oscillations, which required 
constant trim change to control. Flights l and 5 were recorded at 5500 
ft. density altitude in order to avoid turbulence at the 5000 ft. level 
at the tíme of the tests. 

Of interest to individuals involved in this type of flight 
test program is the physical discomfart which may be experienced. The 
authors found that the person acting as recorder is strongly affected by 
nausea, particularly at the short period resonance frequency, while the 
person actually in control of the aircraft is affected but little. This 
leads to the interesting supposition that man is sensitive to third and 


fourth displacement derivatives with respect to time. 





25ο 


DISCUSS ION 


The utilization of frequency response data for an airplane 
has been discussed in the introduction. The primary purpose of this 
investigation was to determine the feasibility of obtaining reliable 
longitudinal frequency response data from steady state dynamic flight 
testing employing simplified instrumentation. With simplicity as a 
primary goal, it was decided to supply a longitudinal forcing function 
to the airplane consisting of a mechanically induced sinusoidal eleva- 
tor deflection of controllable amplitude and frequency. The SFIM re- 
corder provided a lightweight, simplified means of recording this 
elavator input and the airplane's steady state response in normal 
acceleration and pitch rate. Angle of attack response may be obtained 
from normal acceleration response. 

Since the test airplane had a sliding weight installed, it 
was decided to conduct the flight tests at 3 separate C.G. locations 
in order to observe the effect of static stability on the longitudinal 
frequency response. The airplane was weighed in the flight test con- 
figuration while moving the sliding weight through its entire range. 
This enabled a calibration for the airplane C.G. (€ M.A.C.) vs. the 
sliding weight position counter, See Table l and Figs. 7 and B, The 
flight tests were run with the counter readings 0000, 9900, and 979k 
which gave C.G. positions of 25.4, 29.85 and 34.6% M.A.C. respectively. 

In order to observe the effect of flight velocity on fre- 
quency response, it was decided to conduct the flight tests at two 
true air speeds for each C.G. location. The true air speeds used were 
100 and 1,0 M.P.H. The flight tests were conducted at a density alti- 
tude of 5000 feet. 
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In designing the instrumentation system, it was necessary 
to predict the results which could be expected. This enabled an esti- 
mation of the magnitudes and frequencies of elevator oscillation re- 
quired and the approximate mechanical power necessary to drive the 
system. The estimated values of expected pitch rate, normal accelera- 
tion, and elevator deflection were also necessary so that required 
capabilities of the recording components could be determined. 

In calculating the longitudinal frequency ar it was 
necessary to determine the range of oscillation frequencies required, 
Examination of previous tests of this nature indicated that these 
frequencies should include the span from approximately 12 to 150 
cycles per minute. For the flight test configurations to be used, 1t 
was estimated that the period of the phugoid mode would vary fram 
approximately 20 to 30 seconds, Since the period of the minimum ex- 
citation frequency was 5 seconds, it was assumed that little forcing 
of the phugoid mode would be encountered, 

The assumption of negligible excitation of the phugoid mode 
permitted consideration of the steady state flight tests at constant 
velocity. The linearized, longitudinal equations of motion could 
therefore be simplified by eliminating the entire drag equation and 
velocity terms from the lift and pitching moment equations. The latter 
two equations give a simplified solution for the response of the short 
period mode, 

Using these equations, the frequency response for normal 
acceleration ahd pitch rate was calculated as outlined in the calcu- 


lations section of this report. The flight test condition of C.G. at 
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25.4% M.A.Co and V = 1140 m.p.h. was selected as representative of an 
average level of response for the various flight conditions, The sta- 
bility derivatives for the airplane were calculated and compared with 
those obtained experimentally in previous reports, The values consid- 
ered most reliable were used in the calculations, and their values and 
source are shown in the calculation section. The values for stick 
fixed neutral point and elevator power were taken from Princeton Aero- 
nautical Engineering Report No. 275, and are functions of lift coeffi- 
cient. They are plotted in this report as Fig. 9. If it were desired 
to calculate the frequency response for other flight conditions, Xcogo 
could be obtained from Fig. 8, and N, and Un; from Fig. Y. The same 
procedure used in the example solution could then be followed assuming 
that the other derivatives are invariant with flight condition. 

The instrumentation system previously discussed was designed 
using the results of the calculations. A comparison of the calculated 
and experimentally obtained frequency response data for the particular 
flight condition described above is given in Figs. 2l and 25, Although 
there was a sizeable discrepancy, it was extremely gratifying to the 
authors that the original instrumentation design proved completely ade- 
quate throughout the entire flight test range. The capabilities of the 
recording elements were properly utilized and their limits were never 
exceeded. As previously mentioned, a faster film speed for the SF IM 
recorder was highly desirable in order to obtain more accurate phase 
angle determination, but this was beyond the authors' control. The 
deficiencies encountered in the elevator oscillating system were 


quickly remedied. It is felt that this simplified instrumentation 
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system will yield reliable longitudinal frequency response data and is 
relatively trouble free. | 

Samples of flight test data are contained in Figs. 10 and 11. 
The method used for data reduction may be explained by reference to the 
samples of good data contained in Fig. 10. For a particular frequency 
of elevator oscillation, the frequency was obtained by taking 5 cycles 
of the purest sinusoidal traces. The dashed lines at the bottom of the 
film are from the SFIM recorder timer and represent l second intervals, 
With dividers the length of 5 cycles was laid off on the time scale 
such that the time in seconds for 5 cycles was measured, Angular fre- 
quency in radians/sec. was then obtained by mltiplying 5 by 2 Tr and 
dividing by the time in seconds. The amplitude ratios were obtained by 
suming the amplitudes over 5 cycles of each trace with dividers. The 
sumations were measured with a 60 graduations/inch scale. The summa- 
tion for normal acceleration and pitch rate was then divided by the 
summation for elevator deflection. For phase angle measurement the 
length of five cycles was measured on the scale. Then the phase angles 
between the normal acceleration and elevator traces and between the 
pitch rate and elevator traces were summed over 5 cycles and measured 
on the Scale, The particular phase angle was then the ratio of phase 
angle sum to 5 cycles sum multiplied by 3609, The inadequacy of phase 
angle measurement at the higher frequencies is readily apparent from 
observing the distance represented by one cycle or 360? of trace. 

Some examples of poor data are contained in Fig. 1l. In 
each case of poor data, it was readily apparent to the pilots conduct- 


ing the tests and noted on the flight data sheet. In film A, rough air 
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was encountered. As mentioned in the flight test section of this re- 
port, on the flight encountering rough air, further tests were conduct- 
ed at 5,500 feet density altitude. In film B, an example of the pres- 
ence of the phugoid mode may be seen. It was observed during the 
flight tests that occasional gusts and appreciable trim changes would 
excite the phugoid mode, This was evidenced by slow deviations in air 
speed and altitude, and whenever encountered, further flight tests 

were conducted. It is interesting to note that the normal air speed 
variation during the tests was of the order of + 2 m.p.h. which would 
tend to substantiate the constant velocity assumption. In film C, an 
example of erratic pitch rate response may be observed. This resulted 
from an experiment with the mounting support of the rate of turn gyro 
transmittero In this case the gyro was directly mounted on the boiler 
plate without its usual hard rubber pad. In film D, an example of 
elevator distortion may be seen, This resulted from overloading the 
motor by using excessively large elevator displacements, Elevator dis- 
tortion was also present at the lower frequencies before the elevator 
driving motor was geared down. This was due to the motor being over- 
loaded for its driving voltage and the results were an obviously non- 
sinusoidal elevator input. 
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RESULTS 


The frequency response curves of on and = were obtained 
for six flight conditions and are shown on a 12 through 23. It 
should be noted that these plots have abscissas in real time, radians/ 
seco, rather than aerodynamic time, radians/t/_ . The ordinates of the 
pitch rate curves are also in real rather than aerodynamic time. In 
addition, the ordinates of the normal acceleration curve must be cor- 
rected to radian measure if stability derivatives are to be investigated. 

These frequency response curves are a graphical representation 
of the normal acceleration and pitch rate transfer functions of the Ryan 
Navion. Figures 2h and 25 compare experimental and theoretical fre- 
quency response curves for the 10 MPH, 25.h$ m.e.c. flight condition. 
To better understand these results, an analytical expression for the 
experimental results was found by a curve fitting process. The air- 
plane was assumed to have a linear second order characteristic equation, 
D? bD +k=0. The normal acceleration transfer function, neglecting 


tail lift, is then of the form —À o These constants are 
S* + b5 + k 


easily evaluated from plots of amplitude and phase response of a second 
order sinusoidally forced system, 

T ΄ 2 

ου T T E a A 


These plots may be found in any dynamics textbook. The method employed 
to evaluate these constants was an iterative process of selecting a 
trial w,, , finding the gain constant C for this &U,, , then making a «3 
search until the best fit was obtained. From this S anew Wn 

was selected and the cycle repeated. An excellent curve fit was ob- 
tained for the normal acceleration amplitude ratio from the use of 
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CD * 5.95, tf = 0.8, and C = 0.258. A proper curve fit of ampli- 
tude ratio rather than phase angle was chosen because of previously 
discussed difficulties encountered in reducing phase data, The ana- 


lytic expression for the experimental normal acceleration transfer 


function is then Da = 0.258 in real time, or 
s me RO E ea 
δ 86 4 δι86 4 30.3 
ne - 0.200 


er. time. 
= τ...» * 7.158 + 23. in aerodynamic 


The pitch rate transfer function, neglecting tail lift, is 


then of the form SEAL which may be written 9. - Ms (A + Bs), 
SE +S +k 6 3 


The constants A and B were found by the slope intercept 
Vs 

12/5 
expression far the experimental transfer function is then 





method from a plot of ( ) versus frequency. The analytic 


e e 0.258 (6.0 + 2.08) in real tim, or 19 - ασ (5.29 CA 
6 52 « 8.85 4 30.3 6 se + 7.758 + 2305 
in aerodynamic time. Results of these analytic expressions are com- 
pared to the experimental curves in Figures 27 and 28. With the ex- 
ception of the pitch rate phase angle, the airplane appears to be well 
represented by these second order analytic expressions in the frequency 
range investigated. Higher frequency results indicate there may be 
some higher order effects not accounted for. 

The value of the constants b, K, C, A, and B in the experi- 
mental transfer functions are considerably different than calculated 
theoretical values. In addition, a preliminary analysis of stability 
derivatives as determined from experimental data gives some unrealis- 


tic results. Possible sources of error were then consíderedo 
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Repeatability of data eliminates random error as a source. Faulty 
assumptions in the theoretical analysis appear unlikely since proper 
altitudes, air speed, and steady state conditions were achieved during 
the flight tests. Stability derivatives used in the theoretical anal- 
ysis may be a source of error. The static stability derivative used 

in the analysis was obtained from the one "a" flight condition and may 
be quite different from the required En.) Lo value. Another source 
of error could be in the sinusoidal forcing function since no harmonic 
analysis was made. However, examination of the elevator trace indicates 
any error of this type would be very small, 

The possibility of system error in the instrumentation was 
then considered. An error in pitch rate response was found to have 
been introduced by the rate gyro. The frequency response of the rate 
gyro was not considered in the calibration of this instrument, and 
hence proper corrections to rate gyro response were not made in the re- 
duction of flight test data. It had been assumed that since the rate 
gyro has a natural frequency of 5 c.p.S. no resmance in the response 
of the gyro would be experienced. This was probably a valid assumption 
Since the gyro probably has a standard damping factor of 0.7. While 
this would create no error in amplitude ratio response, a considerable 
additional lag due to the gyro is added to the pitch rate response. 

For example, at a forcing frequency of 10 radians/sec. there is a gyro 
lag error of 27° introduced in the pitch rate phase angle response, 
This error accounts for the discrepancy between the analytic expression 
of the transfer function and the experimental curves as seen in Figure 
27. In addition, re-examination of the rate gyro calibration revealed 
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a slightly non linear amplitude characteristic. As much as 10% error 
in pitch rate amplitude could have been introduced due to this non- 
linearity, Other system errors may exist but, if so, they were not 
discovered, 

A major difference in theoretical response and experimental 
response is seen in the characteristic equation. This is not affected 
by the above rate gyro errors, The aerodynamic spring constant k, 
was determined by Wn from the normal acceleration response and is 
felt to be accurate, This term is k > (Cm, +L cmoe) and 
considerable difference must consequently exist Nm these actual 
stability derivatives and the values used for the derivatives in the 
theoretical analysis, Quantitative isolation of this difference is 
felt to be beyond the scope of this report but appears to be & worthy 


subject of future investigation, 
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CONCLUSIONS AND RECOMMENDATIONS 


The method herein discussed for obtaining longitudinal 
frequency response is entirely feasible in most conventional type air- 
planes. The method does, however, have certain inherent disadvantages: 
the most important are as follows: 

lo; Even though the elevator oscillation system used in 

this investigation was uncomplicated, a certain amount of 
Space and rigging was required. This may make the system 
impractical or impossible in many types of airplanes. 

2. The recording instrumentation must have a high degree 

of accuracyo 

3. Smooth flight conditions must prevail. 

h. The method requires considerable flight time. 

So The tests are quite uncomfortable for the pilots. 

In this investigation the data proved to be very repeatable 
indicating little randam error. However, the lack of agreement between 
experimental and theoretical results indicates an instrumentation sys- 
tem error and/or a discrepancy in the theoretical and actual stability 
derivatives must existo The error in the pitch rate phase angle due to 
rate gyro dynamic characteristics is such a system error, The lack of 
correlation between experimental and theoretical characteristic equa- 
tions strongly indicates that the latter error may also exist. 

For further investigations of this type, it is recommended 
that more careful calibration of the instrumentation be made.  Particu- 
lar emphasis should be on the pitch rate gyro to insure, first, that it 


15 linear in the desired range and, second, that the gyro frequency 
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response characteristics be known, It is also recommended that a 
recorder with a higher film speed be used so that accurate phase 
angles can be obtained. 

It is considered that this system is completely adequate for 
further investigation of frequency response and for frequency response 
methods of obtaining stability derivatives. It is recommended that 
further tests be made to positively account for the discrepancy between 
the theoretical and experimental results presented here, 
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CALCULATIONS 


Theoretical approximation to longitudinal frequency response: 
The forcing function used is a mechanically introduced 


sinusoidal elevator deflection of controllable amplitude and frequency: 
Se = O, Sin cot 


The frequency of elevator oscillation, w , will be varied 
to obtain the amplitude and phase angle frequency response for normal 
acceleration and pitch rate. It is assumed that w will be kept suf- 
ficiently high, such that any excitation of the airplane's phugoid mode 
will be negligible. This permits the assumption that V = constant. 
The standard, linearized, longitudinal equations of motion may then be 
Simplified as follows for considering only the airplane's short period 


longitudinal response: 


Lift Θ CL 一 = ÓN 
eq (Es + D) o DO = 
Pitching Mom. eq» — (Cry + Cn, D)X +(Cmpo - AD)DO = -Crs S 
Normal Acceleration Response: 
- Cs § - | 
2 
- me 9 Cm... 7 HD 
E 5 and Cus «(πα 
(x e D) - | SE 


(mat Om D) (Img F: hD) 


S. [-Cng + C5. (Cm - hb) | 5 


ES -hD^« (- ch + (mp, + (npe)D + (Cn + 党 (moe) 
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Let: bz (- s he n, - Cog) 


K = ( Cnx T = (mpg) 


1 , 
i. xX = - (ms [1 E ¿ez (Cm po HD) | 
9 -hD^-bD -K« 


To determine the amplitude ratio and frequency response, 


consider X and 5 as rotating vectorso. Substitute 


D = (Ww into the equations: ( W in radians/, |e ) 


X -Cn 1- 23 t (sez)i | 


-< 
-- 


s (K+hw*) +(bw)e 


- Cni Y hu Y 


15% + (kr hw) 





on Ad 


tam” 2 
á «Ρως 


For the case where ps < f and EA < 1 


the above eq. reduces to: (For example case this applied) 





ο - - Ems tan” e» tan” = m 
9 


Khu 
» 2 
bo^ (K+ hw) 
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— — 
To canvert the above ratio > to = the lift 


equation on the preceding page can be reduced in the following manner: 


Om Cm 
tx x + (Dx -DO) - 2 =o 
e e e ; 


(DX - DO) = non dimensional centrifugal force = - = Ma. 


Nz = incremental load factor positive up = (bad factor -1). Again 


— 


consider ος and δ as rotating vectors and the equation 


becomes ; 
(ως ος "m Cis M" Cm EOD 
EE 2 s Οι 
νο = (Ge a. Cms 
$ λα $ Z4 
_. _ Ems A -H hw A 
Or Na - BEN . .. . Lan EA Mu án Y Mur lmg 
ΕΝ AS Ge 


The above equation will give the amplitude ratio and phase 


Da cu. t be added 
le for —= . The correction factor ( ) mus a 
ο δ Ω Be 


vectorially; in this case, this is a significant factor and must be 


included. 





LO. 
Pitch Rate Response j 


(Use the same procedure as before) 


$e ds 
. (Cg + Gro 0) - (ng 
(Ss, D) zi 


(Cm, top D) Cre | hD) 
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"Ab" rbD = 
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n 0r ao 
b. Viera, (1 Sesh? | (Se BE) nn 
DB. (e) rl zie) 
no 
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A J É 
E LA an σε 


For this airplane (Fe) K Cix and 
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( zen) </ 
as ο } 
therefore, the equation can be simplified to the following : 
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Example Solution: 

v = 140 mh = 205 ft/seo 
ΧΟοβ. = 00.254% M.A.Co. (Counter at 0000. See Figo 8) 
Ao = 5000 feet 
Cy = 2 Wave, m 2(2825) = 0.356 

e sy? (4.00205) (184.2) (205) 
No = 0.364 (Neutral point $ M.A.C.) (See Fig. 9) 
E «ών. 
一 2.61, 
m = Wave. „ 2825 = 87,8 slugs 

g 32.2 : 

T = m = 87.8 = 1.135 sec 

G Y (.00205) (181 .2 (205) 
pP * esi ' το = 10.8 
n = 3,000 slug ft< (An approximation from Princeton Aero- 

nautical Engineering Report No. 231) 

f > = 3,000 = 31.2 1° 

m 87.8 

2 [Eyy2 2 (34.2) 
h = SL 5 = 

A (2) ποτε ης ^ 009 
CL y = 1.81 (Princeton Aeronautical Engineering Report No. 231) 
Cmo "0.08 7) 

Theoretical valves 

pg — 7 70M) 
Ory - ΟἹ. (Log No) = 14.81 (.25h- 36h) = -0.53 
ms = -=.01785 per degree (See Fig. 9) 
b = 


-Cig σ C 
=P de f me „öl 0515) -0.08 -0.17 = -0.3 
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c επ 
K (Cn, En) = -0.53 - 458(0.17) = -o.9ho 


ec - 
E E (2135) gs = -0。237 


BEN) a 


m WEY) e e E ade: UE 3) 
2 2 5.20 Ju + 0 
'$ Yom 0 e (-.940 9.0515 1^ ) 
zs | a.) 
— es hn C 


5 Y(etee2uf-355 ΕΕ 


= Sf Tela 
ὁ . SIS\SB rw κ jan (552) ~ |e" GFE) 8.3 dy” 
A | ως. δω 333 





These calculated plots are shown in Figures 2l and 25 
and are compared to the test data for this case. 
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TABLE 1 


Measurement of Weight and Calculation of Center of Gravity 


Conditions: 
l. Landing gear retracted, flaps up. 
2ο Full gas load. 
3. Canopy closed. 
ho Two pilots with parachutes in forward seats. 
5。 Additional l2-volt battery in aft compartment. 
6. C.G. shifter weight installed. 


7。 Thesis instrumentation in the airplane. 


Measurements : 


Dist. of C.G. | Dist. of C.G. 
Weight | Weight | Total aft of aft of 
C.G. Indicator | Forward Aft | Weight Jack Pts. of MAC 
Full 
0000 Forward 2705 110.5 2815 11.16 25 hh 


6 




























979) Full aft 





Flight Condition 


0.4. = 25.4% M.A.C. 


V = 140 MPH 
HD = 5,000 Ft 


Tal 


4.27 
3.38 
2.20 
1.18 
1.19 
11.89 
4013 
3027 
2.13 
1.91 
l1. 05 
3.19 
2.06 
1.21 
16.31 
13.62 
11.62 
7085 
5.71 
18.46 
13.91 
12.56 
10.68 
10.05 
8.05 
7023 
1.76 
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TABIE 2 


0.201 
0.216 
0.235 
0.253 
0.260 
0.167 
0.195 
0.226 
0.251, 
0.181 
0.203 
0.225 
0.251 
0.256 
0.045 
0.051 
0.062 
0.097 
0.160 
0.038 
0.052 
0.059 
0.071, 
0.081 
0.100 
0.11; 
0.17} 


Llo 








C.G. 
V 


Ho 


Flight Condition 


25 oS MAC. 
100 ΜΡΗ 
5,000 Ft 


TABLE 3 


l5. 


a (ead se | A | ΕΤΕ 


i 
1 
1 
1 
1΄ 
T 
1l 
1Α 
1Α 
1Α 
5 
5 
5 
5 
E 
6 
6 
6 
6 
6 
6 
6 
6 
6 





1.6) 
1.72 
1.75 
1.74 
1.73 
1.61 
1.45 
1.28 
1.62 
1.73 
1.65 
1955 
1.09 
0.80 
0.65 
0.47 
1,98 
1.50 
1.25 
0.98 
0.7h 
0.63 
0.61 
0.52 
0.40 


l6. 
TABLE |, 


Flight Condition 


C.G. = 29.85% M.A.C. 
V = 140 MPH 
Hp = 5,000 Ft 


ee: 


2 
2 
2 
2 
2 
2 
2 
3 
3 
3 
3 
3 
l 
l 
ly 
y 
L 
5 
5 
5 
5 
6 
6 
6 
6 
6 
6 
6 
6 
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TABLE 5 

Flight Condition 

C.G. = 29.85% M.A.C. 

V = 100 MPH 

Hp = 5,000Ft 

ETT EE 1 δαν | Ae | Bis 

2 4.79 0.0813 1.310 
2 11.30 0.0908 2880 1.772 -229° 
2 3.30 0.1155 «2609 1.900 -2069 
2 2.68 0.1335 „251° 1.900 -2009 
2 2.09 0.1472 -2330 1.750 -1920 
2 1.149 0.1680 -2299 0.81, -2009 
3 li. 6l 0.0855 = 1.80 - 
3 n 0.0861 - 1.81, - 
3 3.06 0.1185 - 2.15 S 
3 2.08 0.1295 - 2.01 - 
3 1.75 0.1415 - 1.87 -191° 
ly 5.03 0.062); - 1.75 -2380 
ly 1.23 0.0905 = 2 .0l -2229 
h 3.34 0.1070 - 2.23 -214,9 
lh 2.34 0.1392 - 2.15 -1960 
lh 1.74 0.1610 - 1.91 -1800 
5 18.10 0.0205 - 0.558 - 
5 15.25 0.0212 - 0.600 - 
5 11.90 0.0231 ~ 0。70h - 
5 9.42 0.0271 - 0.880 - 
5 6.72 0.0160 - 1.61 - 
6 17.35 0.0233 - 0.167 - 
6 15.20 0.0255 - 0.540 - 
6 12.50 0.0275 -360? 0.623 -270° 
6 12.20 0.0306 -350° 0.635 -270° 
6 11.00 0.0306 -345° 0.719 -270° 
6 9, 06 0.0367 -31,0° 0.839 -2629 
6 7.35 0.0128 -301,° 1.100 -21,29 
6 6.15 0.0178 -3059 1.390 -237° 
6 1.28 0.0815 ~295° 1,870 
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TABLE 6 


Flight Condition 


C.G. e 3h.6€ M.A.C. 
V =» 10 MH 
Hy = 5,000 Ft 
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Flight Condition 


0.4. e 34.6% M.A.C, 
V = 100 MPH 
Hp = 5,000 Ft 


1.61 
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